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ABSTRACT 

Context. The disks of gas and dust that form around young stars and can lead to planet formation contain polycyclic aromatic 
hydrocarbons (PAHs) and very small grains (VSGs). 

Aims. In this Paper we analyze the mid-infrared (mid-lR) emission of these very small dust particles in a sample of 12 protoplanetary 
disks. Our goal is twofold: first we want to characterize the properties of these particles in disks and see how they are connected to 
interstellar matter, and second we investigate the possibility that their emission can be used as a probe of the physical conditions and 
evolution of the disk. 

Methods. We define a basis made of three mid-lR template spectra: PAH", PAH+, and VSGs that were derived from the analysis of 
reflection nebulae, and an additional PAH' spectrum that was introduced recently for analysis of the spectra of planetary nebulae. 
Results. From the optimization of the fit of 12 star+disk spectra, using a linear combination of the 4 template spectra, we found that 
an additional small grain component with a broad feature at 8.3 fjm is needed. We find that the fraction of VSG emission in disks 
decreases with increasing stellar temperature. VSGs appear to be destroyed by UV photons at the surface of disks, thus releasing free 
PAH molecules, which are eventually ionized as observed in photodissociation regions. In contrast, we observe that the fraction of 
PAH' increases with increasing star temperature except in the case of B stars where they are absent. We argue that this is compatible 
with the identification of PAH" as large ionized PAHs, most likely emitting in regions of the disk that are close to the star Finally, 
we provide a UV-dependent scheme to explain the evolution of PAHs and VSGs in protoplanetary disks. These results allow us to put 
new constraints on the properties of two sources: IRS 48 and "Gomez's Hamburger" which are poorly characterized. 
Conclusions. Very small dust particles incorporated into protoplanetary disks are processed while exposed to the intense radiation 
field of the central star. The resulting shape of the mid-lR spectrum can reveal this processing and be used as an efficient probe of the 
radiation field i.e. luminosity of central star. 

Key words. Astrochemistry — ISM: dust — ISM: lines and bands — Infrared: ISM — Planetary systems: protoplanetary disks — 
Infrared: stars — Methods: observational 

1 . Introduction disk, together with the emission of stochastically heated macro- 

mole cules and nanograins emitting fr o m the surface layers of the 
It IS widely accepted that low-to-intermediate mass (< SM©) ^^^y. dAcke & van den [20041 THabart et al.1 [2004). This 

young stellar objects (YSOs) are surrounded by a disk of emission is characterized by "aromatic infrared bands", gen- 
gas and dust (see e.g. | vanBoekel etal. | |2004| and references .ttrihntPH to nolvr-vHir- arnmatir- bvHrnrarhnns rPAHs 



ga. a..u UU.L ^.cc e.g. | vdUDocKCi ci di. , I f^^^l 'i"" icic cu^c. erallv attributed to polycyclic aromatic hy drocarbons (PAHs, 
therein). There is growing evidence (see [Setiawan et alj|2008| |Leger & Puge j[T984FlAllamandolIitin[T98^ 



and references therein) that these disks are the cradles of plan- 
etary formation, which is why they are called "protoplane- gg^^^^j ^^^^-^^ ^^^^ investigated the properties of PAHs in 

tary disks". From a spectroscopic point of view, these YSOs di sks. Nearly ha lf of the Herbig Ae/Be (2-8 M^ ) stars observed 

are characterized by a large infrai-ed excess due to the emis- bv lMeeus et al.1 (l200l and lAcke & van d^^TAnc ker (2004) show 
sion of dust present in the disk. The near-to-mid-IR region emission. Conc erning their less massive counterparts, T- 

of theii- spectrum is dommated by the emission of hot dust j^^^ ^^^^^ 2M©). lGeers et all (12006,) tend t o show that PAH 

at thermal equilibrium and present in the inner regions of the enu^sion is more marginally detected, while iBouwman et al.1 

(|2008|) report several sources with a strong emission feature 

Send offprint requests to: ^ g 2 urn. These authors, as well as Sloan et al. (2005!), 

O. Berne, e-mail: olivier.berne@cesr.fr , ^, .\. ■ . rr. nAn f . j • c . • 

show that the mid-IR PAH features undergo signmcant varia- 

This work is based on observations made with the Spitzer ^- £ . »u » n v i • 

„ j,,Tr, I T, tions from one source to the next. However, linking these vari- 

Space Telescope, which IS operated by the Jet Propulsion Laboratory, . , ,.n ■ > > , • > ,• • /• 

California Institute of Technology under a contract with NASA. Based ^.^1°"^ ^o the modification of local physical conditions (in par- 

on observations with ISO, an ESA project with instruments funded by ticular the UV field hardness) or PAH properties (e.g. loniza- 

ESA Member States (especially the PI countries: France, Germany, the tion state, size) remains challen ging (see Pe eters et al. 2002t 

Netherlands and the United Kingdom) and with the participation of van Diedenhoven et al. 2004; Sl oan et al.ll2007i on this subject). 

ISAS and NASA. Recently, .Boersma et ah G2008i) have indeed argued that such 
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modifications were linked to the chemical evolution of PAHs in 
the vicinity of HAe/Be stars. Making this connection is essential 
for two reasons. First, it will cast a new light on the properties of 
carbonaceous particles in the interstellar dust cycle, and second, 
this might help in the understanding of disks evolution leading to 
planet formation such as coagulation and evaporation processes. 

Unfortunately such a detailed analysis is difficult given the 
small spatial scales at which these processes occur One way to 
overcome this issue is by increasing the performance of mid- 
IR telescopes and instruments. In this field, impressive obser- 
vational tours de force have bee n performed recen t ly revealing 
the mid-IR anatomy of disks (Laga ge et al.l 120061 : iGeers et alj 
12007; Dou cet et a l. 2007). Another way of proceeding is by tak- 
ing as much advantage as possible of the spectral information. 
Recently, we have studied the link between the shape of the 
observed mid-IR spectrum and the nature of emitting particles 
in various photodissociatio n regions (PDR s ) where the spatia l 
scale is much larger ( Rapacioli et al J [20051: Berne et al. I I2OO7I) . 
We found that their emission spectra could generally be decom- 
posed into 3 fundamental spectra attributed to a population of 
very small grains (VSGs), neutral PAHs (PAH*'), and ionized 
PAHs (PAH^). We were able to link the presence of each one of 
the emitting populations to the local physical conditions: VSGs 
are present in the densest regions of the clouds and dissociated 
into pah" in the more UV irradiated regions to be eventually 
ionized into PAH^ in the cavity around the illuminating stars. 
Following these studies, we propose to analyze the mid-IR emis- 
sion spectra of 12 YSOs showing strong aromatic emission, us- 
ing the set of template spectra of VSG, PAH**, and PAH^ derived 
from our results obtained on PDRs. 

We first describe the observations and the selected sample of 
objects in Sect. |2l In Sect. [3] we justify the choice of template 
spectra and proceed to the fitting of observed spectra. In Sect. |4] 
we discuss on the nature of the diff'erent emitting populations and 
their link with local physical conditions in disks. In particular 
we explore the general trends in the variations of the molecular 
(PAH) versus very small grains emission in disks as a function of 
the luminosity of the central star. Conversely, we illustrate how 
these trends can be used to constrain the luminosity of sources 
that are ambiguous in spectral type and nature, such as IRS 48 
and Gomez's Hamburger 

2. Observations and data reduction 

The selected sources (Table [1]) for this study were taken from 
the large number of Herbig Ae/Be and T-Tauri stars that show 
PAH emission. The objects in this sample are very different from 
one another: some are evolved and probably forming planets 
(e.g. HP 141569 seelAugereauet aljl999: MouiUetet alj200U 
iBrittain & Rettig|[2002h . while others are young and still prob- 
ably enshrouded in the remnants of their parent cloud (e.g CD 
-42° 11721 see Henning et al. 1998). The sp ectral type of IRS 
48 is not determined well , as discussed in iGeers et al.l (l2007h 
and in Sect. 15.11 Furthermore, the nature of Gomez Hamburger 
(IRAS 1805-32 1 1) is debatable. It was originally classified as a 
post-AGB star dRuiz et alj|1987h . but recent observations would 
instead classify it as a young star surrounded byjts protoplan- 
etary disk seen edge-on (see Buiarr abal et al.l 120081 and Hubble 
Space Telescope press release 2002/19). We have only consid- 
ered objects with prominent features detected at 6.2, 7.7, and 
11.3 yum. All the sources are free of the 9.7 jjm silicate emission 
feature, in order to limit the contamination of this band to the 
7.7, and 1 1.3 jum aromatic bands. It is important to note that se- 
lecting disks that lack silicate emission constitutes a bias in our 
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Fig. 1. Mid-IR spectra of the 10 disks in our sample. All are 
from Spitzer-IRS observations except BD -1-40° 4124 and CD 
-42° 11721, which are ISO-PHOT spectra and only cover the 
mid-IR range up to 11 .5 jum. Intensities of the spectrum of BD 
-1-40° 4124 were divided by 10^ for clarity in the figure. 



sample. The absence of silicate emission in the observed spectra 
can be due to several mechanisms: either the silicate grains are 
too far from the central source and thus not heated to tempera- 
tures that are high enough to allow emission in the mid-IR, or 
these grains are too large to be heated efficiently. However such 
effects are (at least at the first order) not related to the nature and 
evolution of carbonaceous grains in disks, so we argue that this 
bias is minor. 

2.1. Spitzer Observations 

HD 135344, HD 169142, HD 34282, HD 97048, HD 34700, RR 
Tau, HD 141159, LkH^ 215, Gomez's Hamburger, and IRS 48 
were observed with the infrared spectrograph IRS (iHouck et al.l 
2004) onboard Spitzer, in the low resolution (^ = 60 - 127) 
mode. We have only considered the data from the short-low (SL) 
module covering wavelengths from 5 to 14.5 fim. We retrieved 



Berne et al.: PAH and VSG emission in disks 
Table 1. Characteristics of the source in the considered sample 



3 



Object 


log(T 


Sp. Type 


Disk evolution 


Age (Myr) 


Ref 


HD 135344 


3.810 


F4Ve 


Transition 


8 ±4 


1,3,11 


HD 169142 


3.914 


A5Ve 


Transition 


61 


1,4,12 


HD 34282 


3.941 


A3Vne 


Thin 


5 


1,5,13 


HD 97048 


4.000 


A0/B9Ve+sh 


Thick 


4±1 


1,6 


HD 34700 


3.774 


GOIVe 


Thick 


n +6 

-3 


1,7 


RRTau 


3.927 


A4e 


Transition 


3+' 

-'-2.6 


1,8,9 


HD 141569 


3.979 


AOVe 


Thin 


5 


1,5 


BD +40° 4124 


4.291 


B2Ve 


Thick+Env. 


1+2 


1,8,14 


CD -42° 11721 


4.470 


BOIVep 


Thick+Env. 


5+^ 


1,10 


LkHtt 215 


4.114 


B7 


Thick+Env. 


1 


14 


IRS 48 


7 


7 


Thin 


4±1 


2 


Gomez Hamburger 


7 


7 


Thick 


7 


15,16 



1: Acke & van den Ancker (2004), 2: Geers et al. (2007), 3: Doucet et al]ri2006) . 4: 
iRamanetal. (2006), 5: Merm et al. (2004), 6: Lagage et al. (2006*), 7: jTorresI j2004|), 8: 
iNatta et al. (1997), 9: Rostopchina ( 1999), 10: Wang & Loonev ( 2007), 1 1: Brown et al. ( 20Q3), 

12: .Gradv et al, (2007), 13 : Pietu et al.. C2003.) , 14: iHillenbrand et al., 41995.) , 15: iRuiz et'S] 

( I1987I) . 16: lBuiarrabaletaLH2008h 



the basic calibrated (.bed) files processed with the S15 pipeline 
from the Spitzer Science Cent er database. The spectra were built 
using the CUBISM software dSmith et al.ll2007 ) and an aperture 
of 4x2 pixels of 1.8" to extract the spectra. 

2.2. ISO observations 

BD +40° 4124 and CD -42° 11721 were observed between 5.8 
and 1 1 .6 fim with the ISO-PHOT instrument onboard ISO, in 
spectroscopic mode ~ 90). We retrieved the highly pro- 
cessed data from the ISO archive and used it directly. 

2.3. Characteristics of aromatic emission from disks 

The spectrum of each one of the sources from the sample is 
presented in Fig. 1. Each spectrum clearly shows the presence 
of aromatic emission at 6.2, 7.7, and 11.3 um. As described 
by ISloan et a l. ( 2005), the spectra resemble those defined as 
Class "B" by Peeters et al. (2002), having their "7.7" /im fea- 
ture shifted towards 7.9 fim. The coolest objects, HD 135344, 
HD 162149, and HD 34700, have the broadest"7.7" jjm feature, 
and it is shifted towards 8.3 fim. 

3. Analysis of the 6-14 //m spectrum of disks 

3. 1 . Fitting fhe mid-IR spectra 

Recently, 'Jo blin et"an ( l2008l) have applied a fitting procedure to 
a sample of mid-IR spectra of planetary nebulae (PNe). The anal- 
ysis employs a basis of 6 mid-IR (6-14 fim) template spectra 
to fit the continuum subtracted spectra of evolved stars: 3 PDR 
components (VSG, PAH°, and PAH+), very large ionized PAHs 
(PAH^) and two broad features centered at 8.2 and 12.3 jum. The 
PDR components were extracted usi ng bli nd si gnal separation 
metho ds as described in lBerne et alj (l2007l) and iRapacioli et alj 
(l2005h . The templates were then generated by fitting these com- 
ponents with Lorentzians after subtraction of the rising contin- 
uum for VSGs. They include all major features except a plateau 
at A > 12yum for VSGs, which was difficult to reproduce due to 
contamination by the H2 line at 12.3 fim and underlying contin- 
uum. The PAH' spectrum was built empirically to account for 



the mid-IR emission observed in PNe and HII regions but wa s 
inspired by quantum chemistry calculations (iJoblin et al.ll2008h . 
As discussed in this paper, there is some uncertainty about the 
intensity of the 11.3 //m band associated with the PAHMem- 
late s pectrum. The set of template spectra used in I Joblin etaP 
I2008h enables the mid-IR emission arising from PNe and HII 
regions to be successfully interpreted both in the Galaxy and 
Small and Large Magellanic Clouds. We therefore apply it here 
to the star+disk systems. 



In our first attempt to fit the observations, we figured out a 
recurrent problem with the fits obtained for the coolest stars: the 
7.7 jum massif could not be reproduced efficiently as shown in 
Fig. 3 for HD 135344. It is clear from this figure that a broad 
feature at ~ 8.3 yum is needed to account for the observed emis- 
sion. What t ype o f grains can be responsible for this emission? 
iSloan et al.l (f2007b suggest that t his shifted band corr esponds to 
the class C population defined bv lPeeters et al.l (l2002i) and is due 
to aliphatic bonds in emitting mixtures. The fitting strategy we 
have used includes a template spectrum defined by a broad fea- 
ture at 8.2 fim (called "BF" in Joblin et al. 2008), based on a 
class C spectrum and prominent in post-AGB object. The car- 
riers of this feature are observed to be easily destroyed by UV 
photons, therefore their detection in protoplanetary disks would 
imply that they are reformed in the parental molecular cloud or 
in the densest environments of the disk itself. Since the chemi- 
cal conditions differ from those prevailing in evolved stars there 
is no reason why these species should be identical in both en- 
vironments, and this explains why the 8.2 jjm BF is not able to 
reproduce the emission of HD 135344. We therefore define a 
new BF template spectrum and refer to it as dBF (for disk broad 
feature) to avoid confusion with the band seen in evolved stars. 
The dBF template spectrum is based on the observed spectra of 
T-Tauri (Fig. 14 in Bouwman et al. 2008). The new fit of HD 
135344 is presented in Fig. |4] Following J obhn et all (|2008|) we 
also used 2 versions of the PAH^ spectrum. Results do not vary 
significantly whether the 11.3 and 12.7 jjm features are present 
in the template spectrum or not. This implies that the shape of 
the PAH^ spectrum is not well-constrained in the 11.3/12.7 jim 
range. It is not excluded that the carriers of the dBF also have 
some features in the 10-14 yum range, but these are expected to 
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Fig. 2. The five mid-IR template spectra used to fit the observed 
emission of protoplanetary disks. The PAH", PAH^, and VSG 
spectra are PDR componen ts. The PAH'^ c omponent (very large 
ionized PAHs) is from Joblin et al.1 (l2008h and the 8.3 jum dBF 
was introduced in this paper. 



be weak. In the following, we concentrate on the results includ- 
ing the 11.3 and 12.7 fim feature in the PAH^ spectrum, and no 
additional fe ature in the 8.3 fim dBF spectrum. As discussed in 
iJoblin et all ((2008), the main effect of using a PAH* spectrum 
without 1 1.3 and 12.7 /urn feature will be to slightly decrease the 
fraction of PAH° in the results of the fits. Finally, we note that 
the 12.3 yum BF used in lJoblin et alj ((2008) is useless in the fits 
of protoplanetary disks, so we do not use it here. The final basis 
of 5 template spectra we have used for the fitting procedure is 
presented in Fig. [2] 

Considering the geometry and density of disks, IR radiative 
transfer might be a concern. We have therefore re-run our fit us- 
ing simple extinction parame ter based on the extinction curve of 
IWeingartner & Draind (1200 lb . for a total-to-selective extinction 
ratio Rv=5.5. Though in some cases the fit is slightly better, the 
estimated fraction of each population is nearly unchanged. Thus 
, for the sake of simplicity, we present here results that do not 
include extinction. 




X (jim) 

Fig. 3. Preliminary fit of HD 135344, using only PDR com- 
ponents together with PAH'. Observed spectrum is the upper 
dashed line and fit is the continuous line. The contribution of 
VSG, PAH^, and PAH*' is represented below in black, red, and 
green, respectively. It appears clearly that the use of these com- 
ponents cannot reproduce the observed emission and that a new 
component with a feature at 8.3 yum is needed (see dBF in Fig.[2]i. 
The new fit is presented in Fig.[4[ 



3.2. Results 

The results of the fits are presented in Fig. [4| For each object 
we are able to reproduce the observed spectrum quite efficiently. 
Table [2] gives the fractions of mid-IR integrated fluxes for each 
component present in the studied source. Because mid-IR emit- 
ters are heated and processed by the UV-visible photons emerg- 
ing from the central star, we propose to compare the fractions 
of each emitting population to the theoretical total and UV lumi- 
nosities of the central source given its spectral type (see Table 1). 
We used the spectra from the 1993 Kurucz stellar atmosphere at- 
las ( http://www.stsci.edu/hst/observatory/cdbs/k93models.htrni' 
based on lKuruczlll979i) . We used a gravity of 4 for all sources 
and a stellar radius from lAller et al., (.1982) . The total lumi- 
nosity is defined as A-nRt J F^dA and the UV luminosity as 

47tR, ^-''^^ (according to lHabing|[T968[) . where Rt is the 

stellar radius and F,i the flux from the Kurucz catalog. In Fig. [5] 
we plot the fraction of VSG emission as a function of UV lumi- 
nosity, showing a clear decrease of this emission while the mass 
of the star increases. In contrast, the fraction of PAH emission 
increases with the luminosity (Fig. (61(71). This is due to the in- 
creasing PAH' emission, whereas the PDR-type PAH emission 
is found to increase up to spectral type A4 and then to decrease 
when going towards stars of earlier spectral types (Fig.O. The 
same trends are observed either using total or UV luminosity. We 
also observe that there are no dBF in the spectra of disk around 
stars of spectral type earlier than AO (i.e. UV luminosity >~ 10^^ 
W). 

4. Probing disks with mid-IR emission ? 

4.1. The origin of tlie 8.3 /vm broad feature 

The 8.3 fim dBF component is found primarily in the sources 
with the coolest stars, which suggests a fragile material as in 
post-AGB stars, possibly of aliphatic nature. Sloan et al. (200^ 
suggest that the grains responsible for the 8.3 fim feature in 
cool protoplanetary disks ar e similar to the grai ns responsible 
for the"class C" spectrum of iPeeters et al.l (|2002|) and observed 
in post-AGB stars. However, we find that the dBF position and 
width are not compatible with the class C spectrum, the dBF be- 
ing broa der and peaking a t 8.3 rather than 8.2 fim, as also noticed 
by .Bouwman et al.l(l2008h . 
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Fig. 4. Fits of the 12 spectra of the star+disk systems observed with Spitzer-IRS and ISO using our set of 5 template spectra. The 
observed spectrum is the dashed black line, the fit is the continuous black line. The contribution of each component is shown below 
in colors (using the same color code as in Fig. 2.) 



4.2. Processing of VSGs 

The fraction of VSG emission clearly decreases with increas- 
ing UV luminosity (see Fig. |5]l, while emission from free 
PAH molecules becomes more and m ore prominant . This 
is in line with observations of PDRs (ICesarsky et al.l l2000t 



iRapacioU et al ] 120051) where it appears that VSGs are dissoci- 
ated into free PAH molecules. It is expected that such a pro- 
cess is also taking place on the surface of the disk. Still, one 
has to con sider the competition between evaporation and co- 
agulation jRapacioli et al. 1 120061) . but evaporation is likely to 
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Table 2. Factions of each population of the carbonaceous dust 
emission in the spectrum of the 12 sources of the selected sample 



fraction of tlie emission (%) 



Object 


VSG 


PAH" 


PAH+ 


PAH' 


dBF 


HD 135344 


66 





1 


24 


9 


HD 169142 


38 


25 


11 


26 





HD 34282 


30 


13 





51 


6 


HD 97048 


37 





10 


50 


3 


HD 34700 


61 


14 





18 


7 


RRTau 


27 


37 





36 





HD 141569 


34 





12 


51 


3 


BD+40 4124 





61 


29 


10 





CD-42 11721 


8 


57 


23 


12 





IRS 48 


27 





13 


58 


2 


GoHam 


54 





32 


14 





LkHQ'215 


25 


43 


20 


12 






dominate when the central star becomes hotter and therefore 
produces more and more UV photons. Dullemond et al. (2007) 
show that in more evolved disks the ratio of PAH-emission over 
far-infrared emission should be much higher than what is ob- 
served, because of the lower optical thickness of the disk in the 
wavelength range at which PAHs absorb. Thus they propose that 
coagulation of PAHs could be the reason for this lack of emis- 
sion. Rapacioli et al. (2005, 2006) suggest that VSGs are indeed 
coagulated PAHs. In this scenario, an increase in the abundance 
of coagulated PAHs would imply an increase in VSG emission 
in more evolved disks which is not observed. Instead, one could 
claim that this lack of PAH emission is due to the geometrical 
evolution of the disks. Indeed, an older and flatter disk will re- 
ceive less UV photons per unit of surface, and thus the heating 
of PAHs (if still present) might become less efficient leading to 
a less intense mid-IR signal. 

4.3. Evolution of PAH populations 

Figure |8] presents the fraction of PDR-type PAHs as a function 
of the stellar UV luminosity. There is a clear rise in the frac- 
tion of emission from sources of spectral type GO to A4, and 
then this fraction decreases again with increasing luminosity. 
The rise can be interpreted as due to the high production rate 
of PAHs from VSGs (see previous section) in the disks around 
stars of GO to A4 spectral types. For sources that are hotter, these 
molecules start being destroyed by the strong UV-visible radia- 
tion field, which explains the decrease in the observed PDR-type 
PAH emission fraction. In contrast, the fra ction of PAH ' emis- 
sion increases (Fig. |6]l. iJoblin et"aLl (l2008h conclude that this is 
the likely result of the destruction of the smallest PAHs com- 
bined with efficient heating of the PAH '' population. This is con- 
sistent with PAH^ being very large ionized PAHs that can better 
survive in extreme irradiation conditions. 

The charge state of PAHs in disks is determined by the com- 
petition between ionization rate versus electron recombination 
for cations/neutrals and photodetachment versus attachment rate 
for anions/neutrals. In disks, the gas density can be as high 
as 10^"^ cm"^, and the electron attachment to PAHs can be 
more efficient than re-ionization by UV radiation as shown by 
[Li & Lunine (2003). This would imply that PAHs can be found 
as anions (PAH ) that could contribute strongly to the mid-IR 
emission spectrum of disks. In this case it is tempting to identify 
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Fig. 5. Fraction of VSG emission in the selected sources, as a 
function of the UV luminosity of the central star Diamonds are 
for sources that have a well-defined spectral type. The spectral 
types of GoHam and IRS 48 (shown with star symbols in this 
plot) are inferred in order to agree with the found correlation. 
The correlation shows that VSGs are less abundant around hot 
stars, likely due to their rapid destruction in these harsh environ- 
ments. 
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Fig. 6. Fraction of PAH (PAH^H-PAH+H-PAff) emission in the 
selected sources, as a function of UV luminosity. Diamonds are 
for sources that have a well-defined spectral type. The spectral 
types of GoHam and IRS 48 (shown with star symbols on this 
plot) are inferred to agree with the found correlation. 



PAH' as very large PAH . However IVisser et al.l (l2007l) show 
that in flared optically thick disk (e.g. HD 97048), and even 
though anionic PAHs could be abundant, they are likely to be 
situated in deep regions of the disk where they receive too little 
UV light from the star to be excited so they dominate the mid-IR 
emission. On the other hand, they suggest that the emission at 
the surface of HAe/Be disks will arise from a mixture of very 
large (> 100 C atoms), positively ionized PAHs. 

Figure|7]shows that the fraction of PAH' emission increases 
with UV luminosity for the studied T-Tauri and Herbig Ae stars. 
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For the Be star nearly no PAH' emission is detected, which 
we interpret as due to the absence of disk emission in the ob- 
served spectra (see Sect. 14.5b . so it is not plotted on Fig. |7] 
Furthermore, these objects do not exhibit the dBF and their spec- 
trum can be reproduced using the PDR spectra almost exclu- 
sively. Spitzer-IRAC images as well as millimeter observations 
(lHenningetal.lll998h suggest that these sources are not isolated, 
and the ir mid-IR spectrum wa s indeed previously classified as 
"A" by Ivan Diede nhoven et alj (|2004), which is consistent with 
typical PDRs. This suggests that we are in fact not observing the 
emission from the disk but fr om the PDR created in the nebulos- 
ity around the central source. lBoersma et al. 1(12008) come to the 
same conclusions concerning the B9 stars HD 37411 and HD 
36917 observed with Spitzer, for which they suggest that the 
observed emission is arising from the surroundings of the star 
but not a disk. There can be three reason for this: (1) no PAH 
can resist in the disk because they are all photodissociated, (2) 
there is no disk because it has been already photoevaporated, (3) 
the geometry of the disk is such that the PAH emission is not 
sufficiently activated by UV-visible photons. We discuss these 
scenarios in Sect. 14.41 

That some PAH' emission is present in cool stars, as in the 
F4 star HD 135344, implies that the near-UV radiation field 
is sufficient to destroy some of the smallest PAHs and induce 
PAH"^ emission. This means that PAH' are close to the star de- 
spite HD 135344 havin g a rather large inner radius (see e.g. 
iPontoppidan et al.ll2008l) . In this case, one could consider that 
PAHs are in fact coupled w ith a gas disk that extends closer to 
the star than the thick disk (IPontoppidan et al] [2008). The pres- 
ence of PAH' in HD 135344 argues against the possibility that 
these species are positively ionized as the F4 star is unlikely to 
produce enough UV to ionize the PAHs. On the other hand, PAH 
cations appear to be good candidates considering that they are 
very abundant in emission in the disk aro und A star s that p ro- 
duce a large amount of UV (Fig.|2ll. As in Joblin et al. (2008) we 
conclude that PAH^ are very large ionized PAHs tracing harsh ir- 
radiation conditions. Unfortunately, the observational constrains 
provided here are not sufficient to determine if PAH ' correspond 
to a population of cations or anions or even a mixture of both. 
We emphasize nevertheless that the detection of PAH ' around a 
young star provides a way to assess the presence of a disk. 

4.4. The mid-IR spectrum: a probe of the irradiation 
conditions of the disk 

Finally, can we trace the morphology/content of the disks using 
the spatially non-resolved mid-IR Spitzer spectra? The answer 
is probably no, because the evolution of very small dust parti- 
cles in disks is clearly connected to the irradiation conditions 
and thus to the spectral type of the central star. This connection 
between the composition of the mid-IR spectrum and spectral 
type is presented in a schematic way in Fig. |9] and corresponds 
to the following scenario: the dense disk is a reservoir of VSGs 
(such as molecular clouds in PDRs), which are constantly disso- 
ciated into free PAH molecules. The more intense the UV field 
(i.e. the more massive the star is), the shortest the lifetime of a 
VSG at the surface of the disk, where the emission is occurring. 
This is why the fraction of VSG emission decreases with inreas- 
ing star temperature (Fig. |9]l- Similarly, once they are released 
from VSGs, the lifetime of the smallest PAHs around hot stars 
is too short to allow their emission to be dominant. This is why 
the fraction of large PAHs (PAH^) that resist destruction more 
is greater in disks around hotter stars (Fig. |9]l. More massive B 
stars only show emission of PDR-type PAHs. This might be be- 
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cause this emission does not arise predominantly from a disk, 
either because VSGs and PAHs (even very large PAH') do not 
survive in disks around Be stars or because the disk is small or 
absent (see next section). 

4.5. Photoevaporation of disks around Be stars 

We have shown that Be stars have a PDR-type emission and that 
the PAH^ emission is not observed in these objects. This sup- 
ports the idea that there are no disks around these stars or that 
the disk emission in the mi d-IR is too weak to be disentangled 
from the surrounding PDR. IVisser et al.l (l2007h have shown (see 
their Fig. 3) that a radiation field of more than ~ 5 • 10^ in units 
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Fig. 9. Schematic representation of the contributions of the very 
small dust particle populations to the mid-IR spectra of T-Tauri 
and Herbig Ae/Be stars. While the spectrum of G, F, and A stars 
is dominated by disk emission, B star spectra are likely due to 
emission arising from the remnants of the parental cloud. The 
total emission is due to the contribution of VSGs (dashed line), 
PDR-type PAHs (dash-dotted line), and PAH' (dotted line). 



5. Putting new constraints on the nature of IRS 48 
and Gomez's l-iamburger 

5.1. The luminosity of IRS 48 

The spec tral type of IRS 48 is controversial (see discussion in 
iGeers et al. 2007). Here we propose to classify it using the re- 
sults of our analysis. With the observed abundance of PAH*, 
VSGs, and the presence of the dBF, we can estimate in which 
range of temperature IRS 48 falls. The PAH' fraction in the 
emission of IRS 48 is 58%, which is compatible with a luminos- 
ity of ~ 5 • 10^^ W based on Fig. 4. The VSG fraction is 29 %, 
which also compatible with a UV luminosity of ~ 5 ■ 10^^ W (Fig. 
5). Finally, the fraction of dBF in the spectrum of IRS 48 is 2%, 
which puts it near the limit for the presence of the dBF, consis- 
tent with a UV luminosity of ~ 10^^"^^ W. Further more, we ob- 
serve a 7.7 jum to 1 1 .3 fim band ratio above 2. Using .Flagev et al.] 
(l2006h . this implies a — above 10^. Using a temperature T of 

500 K a nd an electron density of ne=10'^ S."^^ surface of 

the disk (iJonkheid et al.ll2007t IVisser et al ] l2007h . this leads to a 
UV field of Go ~ 60000 in units of the Habing field, which is 
similar to what is found in the emitting regions of disks around 
an AO star. This value, as well as the derived luminosities, are 
compatible with the powering source being an early A star. If 
the spectral type of the source is really around F3 or even MO 
(Geers et al. 2007), an excess of UV is needed and might indi- 
cate accretion processes. 



of the Habing field will destroy even large PAHs of 100 carbon 
atoms in less than 3 Myrs. An estimate of the radiation field at 
a distance of 300 AU from a B2-0 star, taking only dilution ef- 
fects into account and no UV radiative transfer, gives values as 
high as 10^"^. This is high above the limit value of 5 ■ 10^, so it 
is likely that in such a harsh environment most, if not all, PAHs 
will be destroyed efliciently at a radius of 300 AU and even more 
efficiently at smaller radii. A second possibility is that the disk 
itself is absent because it was photoevaporated. It is indeed be- 
lieved that UV can provoke the dispersion of the gas and small- 
est dust particles (therefore probably including PAHs and VSGs) 
within a few 10^ Myr (see e.g. the r ecent review by Alexander, 
2007). Recently, lAlonso-Albi et al.l (|2008) have found that the 
dusty disks around Herbig Be stars tend to be lighter and smaller 
than those towards cooler stars, probably due to the photoevap- 
oration efi'ect. The PAH emission in such a small disk is thus 
expected to be very weak in any case (whether or not PAHs are 
destroyed). Note that the eff'ect of PAH destruction in the disks 
around Be stars will also directly aff'ect photoevaporation of the 
disk. PAHs and VSGs are responsibl e for most of the UV-optical 
thickness of the disk ( lDrainell2003h so their destruction would 
imply a decrease in the UV optical thickness. This effect com- 
bined with the higher proportion of UV photons, will tend to 
accelerate the photoevaporation process of the disk and implies 
that the lifetime of disks around massive stars could be shorter 
than the ti me needed to dispe r se the surrounding material as sug- 
gested bv lAlonso-AIbi etakl (l2008l) . The probably flatter shape 
of the disk around Herbig Be stars (because of the photoevapo- 
ration of the outer layers) could also contribute to the absence of 
PAH emission because of the self-shadowing eff'ect. 



5.2. The nature of Gomez Hamburger 

Gomez Hamburger (GoHam hereafter) is subject of debate in 
term s of classification . It was originally classified as a post-AGB 
star dRuiz et al.lll987h . and recent observations would instead 
classify it as a pre-main-seque nce (PMS) A star surrounded by 
its disk (Buiarraba l et alj|2008l) . We propose to investigate here 
the observed mid-IR spectrum of GoHam is compatible with the 
PMS-star + disk hypothesis. As described above, we can use the 
results of the fitted spectra to investigate the properties of the il- 
luminating star The found fractions of VSGs and PAH ' and dBF 
carriers are compatible with a total luminosity of ~ 8 Lq and a 
UV luminosity of ~ 0.3 Lq. This indicates that the central star 
should have an effective temperature of ~ 7000-9000 K. This is 
slightly cooler than the stellar temperature of ~ 10 000 K de- 
rived by Ruiz et al. (1987) from near-IR photometry. However, 
one should note that measuring the temperature of the star us- 
ing photometric bands is extremely difficult and unreliable given 
that the disk is seen edge on, and thus that the observed light is 
completely scattered by the disk. Based on our analysis, the cen- 
tral star in GoHam would be an A5-8 star near main sequence 
(n o emission fines a re present in the visible spectrum observed 
bv' Ruiz et al.lll987t) . This imp fies a stellar mass of about 1.7-2 
Mq from evolutionary tracks dvan den Ancker et al.lll998l) that 
is compatible with the mass derived from CO observations. Our 
derived luminosity of ~ 8 Lq implies a distance of ~ 200-300 pc 
to reproduce the SED. T he Hubble Space Teles cope images, as 
well as CO observations (Buiarr abal et al.ll2008l) . suggest a disk 
radius of 4". Using the distance we derive of 200 pc this gives a 
physical radius of 800 AU, consistent with current estimates of 
the ou ter radius for a disk around HAeBe objects (see compila- 
tion in lHabart et al.l (|2004|) and references therein). Furthermore, 
the A5-8 scenario does not necessitate the presence of a binary 
as suggested by Buiarrabal et al. (200^. 

Finally we note that the fit do es not require the 8.2 //m BF 
commonly seen in post-AGB stars jJoblin et al.l2008l) and shows 
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some PAH^ emission. Thus, if not a PMS-star, GoHam would 
instead be a planetary nebula that can lack the 8.2 jum feature 
and show the 7.9 fiia PAH* feature in their mid-IR spectrum. 

6. Summary and conclusion 

We have studied the mid-IR emission spectrum of 12 disks 
around young stars, showing clear signatures of stochastically 
heated grains and no silicate emission. The observed spectra, 
emanating from very different objects, can be reproduced using 
a mixture of only 5 template spectra, 3 of which are represen- 
tative of interstellar (PDR) populations. In addition, a PAH^ 
population characteristic of highly UV-irradiated environments 
is present as is a peculiar feature (dBF) possibly due to aliphatic 
material in analogy with the 8.3 //m broad feature present in 
post-AGB stars. We can recapitulate the main results of this 
work in three points: 

(1) The mid-IR spectrum of disks traces the PDR created 
at the surface of the disk. As in PDRs, we find that VSGs are 
destroyed to become free PAHs. VSGs and small PAHs survive 
long enough to be observed only in the disks around cool stars. 
Around A stars only large PAHs have a Ufetime that is long 
enough to allow their detection. The coagulation of PAHs at the 
surface of disks, i.e. mid-IR emitting region, is likely hindered 
by the action of the UV field and therefore not correlated with 
the coagulation of larger grains. 

(2) PAH* emission indicates the presence of a disk. This 
emission is found to be more prominant in HAe stars with high 
UV fields. This is compatible with PAH* being very large PAHs 
that can resist photodestruction better. Spectroscopy indicates 
that they are charged and models favor cations rather than 
anions. The PAH* and dBF carriers are absent around Be stars, 
which we interpret as the absence of emission coming from a 
disk. This is Ukely because the disk is small or event absent as it 
was photoevaporated. 

(3) The shape of the mid-IR PAH spectrum of Herbig Ae/Be 
and T-Tauri stars is not related to disk morphology but to UV- 
irradiation conditions. This connection can be used to infer the 
spectral type / mass of the central source. Using this result, we 
show that IRS 48 is of early A spectral type or produces a large 
amount of UV photons in excess due to accretion. Similarly, we 
constrain the nature of the central source of GoHam and show 
that the mid-IR emission spectrum of this source is compatible 
with a protoplanetary disk with a central star of spectral type 
around A5-8. 

Finally, we emphasize that detailed studies of a large sample 
of disks with PAH emission are needed to provide deeper in- 
sights into the evolution of these very small dust particles along 
with the evolution of the disk. Furthermore, spatially resolved 
spectro-imagery, in the aromatic emission range, of such objects 
will be crucial for the understanding of local evolution of the 
disk, leading to planet formation. 

Acknowledgements. This work was supported by the French National Program, 
Physique et Chimie du Milieu Interstellaire which is gratefully acknowledged. 

References 

Acke, B. & van den Ancker, M. E. 2004, A&A, 426, 151 

AUamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, ApJ, 290, L25 



Aller, L. H., Appenzeller, I., Baschek, B., et al. 1982, Landolt-Bornstein: 
Numerical Data and Functional Relationships in Science and Technology 
- New Series " Gruppe/Group 6 Astronomy and Astrophysics " Volume 2 
Schaifers/Voigt: Astronomy and Astrophysics / Astronomic und Astrophysik 
" Stars and Star Clusters / Sterne und Sternhaufen (Landolt-Bornstein: 
Numerical Data and Functional Relationships in Science and Technology) 
Alonso-Albi, T, Fuente, A., Bachiller, R., et al. 2008, ArXiv e-prints 
Augereau, J. C., Lagrange, A. M., Mouillet, D., & Menard, F. 1999, A&A, 350, 
L51 

Berne, O., Joblin, C, Deville, Y, et al. 2007, A&A, 469, 575 

Boersma, C, Bouwman, J., Lahuis, R, et al. 2008, ArXiv e-prints, 804 

Bouwman, J., Henning, T., Hillenbrand, L. A., et al. 2008, ArXiv e-prints, 802 

Brittain, S. D. & Rettig, T. W. 2002, Nature, 418, 57 

Brown, J. M., Blake, G. A., Dullemond, C. R, et al. 2007, ApJ, 664, L107 

Bujarrabal, V., Young, K., & Fong, D. 2008, ArXiv e-prints, 803 

Cesarsky, D., Lequeux, J., Ryter, C, & Gerin, M. 2000, A&A, 354, L87 

Doucet, C, Habart, E., Pantin, E., et al. 2007, A&A, 470, 625 

Doucet, C, Pantin, E., Lagage, R O., & Dullemond, C. R 2006, A&A, 460, 117 

Draine, B. T. 2003, ARA&A, 41, 241 

Dullemond, C. R, Henning, T, Visser, R., et al. 2007, A&A, 473, 457 
Flagey, N., Boulanger, F, Verstraete, L., et al. 2006, A&A, 453, 969 
Geers, V. C, Augereau, J.-C, Pontoppidan, K. M., et al. 2006, A&A, 459, 545 
Geers, V. C, Pontoppidan, K. M., van Dishoeck, E. E, et al. 2007, A&A, 469, 
L35 

Grady, C. A., Schneider, G., Hamaguchi, K., et al. 2007, ApJ, 665, 1391 
Habart, E., Natta, A., & Kriigel, E. 2004, A&A, 427, 179 
Habing, H. J. 1968, Bull. Astron. Inst. Netherlands, 19, 421 
Henning, T, Burkert, A., Launhardt, R., Leinert, C., & Stecklum, B. 1998, A&A, 
336, 565 

Hillenbrand, L. A., Meyer, M. R., Strom, S. E., & Skrutskie, M. F 1995, AJ, 
109, 280 

Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 
Joblin, C., Szczerba, R., Berne, O., & Szyszka, C. 2008, A&A, 490, 189 
Jonkheid, B., Dullemond, C. P, Hogerheijde, M. R., & van Dishoeck, E. F. 2007, 

A&A, 463, 203 
Kurucz, R. L. 1979, ApJS, 40, 1 

Lagage, R-O., Doucet, C, Pantin, E., et al. 2006, Science, 314, 621 
Leger, A. & Puget, J. L. 1984, A&A, 137, L5 
Li, A. & Lunine, J. 1. 2003, ApJ, 594, 987 

Meeus, G., Waters, L. B. F M., Bouwman, J., et al. 2001, A&A, 365, 476 
Men'n, B., Montesinos, B., Eiroa, C., et al. 2004, A&A, 419, 301 
Mouillet, D., Lagrange, A. M., Augereau, J. C, & Menard, F. 2001, A&A, 372, 
L61 

Natta, A., Grinin, V. R, Mannings, V., & Ungerechts, H. 1997, ApJ, 491, 885 
Peeters, E., Hony, S., Van Kerckhoven, C., et al. 2002, A&A, 390, 1089 
Pietu, v., Dutrey, A., & Kahane, C. 2003, A&A, 398, 565 
Pontoppidan, K. M., Blake, G. A., van Dishoeck, E. P., et al. 2008, ArXiv e- 
prints, 805 

Raman, A., Lisanti, M., Wilner, D. J., Qi, C, & Hogerheijde, M. 2006, AJ, 131, 

2290 

Rapacioli, M., Calvo, F, Joblin, C., et al. 2006, A&A, 460, 519 
Rapacioli, M., Joblin, C., & Boissel, P 2005, A&A, 429, 193 
Rostopchina, A. N. 1999, Astronomy Reports, 43, 113 
Ruiz, M. T, Blanco, V., Maza, J., et al. 1987, ApJ, 316, L21 
Setiawan, J., Henning, T, Launhardt, R., et al. 2008, Nature, 451, 38 
Sloan, G. C., Jura, M., Duley, W. W., et al. 2007, ApJ, 664, 1 144 
Sloan, G. C., Keller, L. D., Forrest, W. J., et al. 2005, ApJ, 632, 956 
Smith, J. D. T., Armus, L., Dale, D. A., et al. 2007, PASP, 1 19, 1 133 
Torres, G. 2004, AJ, 127, 1187 

van Boekel, R., Waters, L. B. F. M., Dominik, C, et al. 2004, A&A, 418, 177 
van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 1998, A&A, 330, 

145 

van Diedenhoven, B., Peeters, E., Van Kerckhoven, C, et al. 2004, ApJ, 611, 928 
Visser, R., Geers, V. C., Dullemond, C. P, et al. 2007, A&A, 466, 229 
Wang, S. & Looney, L. W. 2007, ApJ, 659, 1360 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 



